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Abstract 

Einstein equations are addressed with the energy-momentum tensor that appears 
if the equations under discussion are required to possess conformal invariance. It is 

! proved that thus derived equations (equations of conformally invariant geometrody- 

namics) can have not only smooth solutions, but also solutions with discontinuities 
on space-like hypersurfaces. The solutions obtained are similar to the well-known 

1 ^ , discontinuous Einstein equation solutions like shock-wave solutions, extended-source 
solutions, etc. 

■ For the centrally symmetric stationary solution discussed in the paper, the dis- 
continuity surface removes the singularity. The degree of generality of this solution 
regularization mechanism is discussed. 

f--. The issue of the mechanism that forces any smooth solution in the conformally 

invariant geometrodynamics to be rearranged into the discontinuous one when cer- 
tain conditions are met is also discussed. The conditions can be: (1) sound speed 
becoming to be higher than light speed; (2) the solution becoming intolerant to 

■ smaller and smaller-scale perturbation modes. 
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This paper is a summary of ref. [1] that contains missing proofs of the statements 



presented here and needed references. 



We consider discontinuous solutions to Einstein equations with the nonzero energy- 
momentum tensor. The equations are 

Ra/3 — \9apR = T a p, (1) 

where 

T a p = -2A a Af3 - g a pA 2 - 2g aP A v v + A a , p + A p . a + g af} \, (2) 

A a is a so-called gauge vector, A is, generally speaking, one of the desired functions. 
Equations (1) retain their form in conformal transformations 

SW3 -> 9ap ■ e 2fT , A a -> A a - a. a , A -> A • e~ 2a . (3) 
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Therefore we refer to equations (1) with energy- momentum tensor (2) as equations of 
conformally invariant geometrodynamics. 

The Riemannian space complemented with A a , A and transformation rules (3)can be 
considered as Weyl space. In the Weyl space the connectivity is expressed in terms 
of Christoffel symbols (^^j and vector A a with relation 

= (ap) + ^Ap + 5 x p A a - g a( ,g XT A T . (4) 

Equations (1) with energy-momentum tensor (2) possess a unique property: Cauchy's 
problem for these equations is posed without any connections to Cauchy's data. Thus, 
in the synchronous frame (g ok = 0) with using gauge condition A = Const the Cauchy's 
data are the following sixteen functions: g mn , g mn , goo? A k . 

The starting axiomatics of the Weyl space requires that the minimum necessary geo- 
metric object smoothness classes specified in Table 1 be ensured. 

Table 1 - Minimum smoothness classes 



Object 


Smoothness class as 


Smoothness class on 




a whole on the 


local charts (pieces) 




manifold 




Components of metrics g a p 


C 1 


C 2 


Christoffels 


c° 


C 1 


Components of gauge 


Piecewise 


c° 


vectorA a 


continuous 




Function A 


Piecewise 
continuous 


c° 



In this paper our interest is not with light-like discontinuity surfaces, but with space- 
like hypersurfaces. Such discontinuity surfaces are known to appear only when the energy- 
momentum tensor is nonzero. These solutions include gas-dynamic shock-wave type so- 
lutions, extended-source solutions, etc. What follows will be analogous to the shock wave 
theory within the framework of the general relativity theory, but under the sole condition: 
the energy-momentum tensor will be of form (2). 

Apply the requirement of the existence of only minimum smoothness classes to the 
static centrally symmetric problem. The squared integral for this solution, without loss 
of generality, can be reduced to 

ds 2 = -exp (7) • dt 2 + exp (-7) • dz 2 + exp ((3) ■ {d6 2 + sin 2 6 ■ dy? 2 ) (5) 
and the gauge vector to 

A, = (0,0, 0,0). (6) 

For four desired functions 7, /3, <f), A relations (1) result in four equations which we 
will not present here. Instead, we at once present the general solution branch, which can 
transfer to de Sitter solution and Schwartzschield solution. 

<f) = p ■ exp (7) 
exp {(5) = A ■ sh 2 (pz + a) 

exp (7) = -^j + B ■ [pz ■ cth (pz + a) - 1] + bp ■ cth (pz + a) 
A (z) = Bp 2 



(7) 
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Seek the solution, such that has a discontinuity for radial variable z — A. In region 
z < A the desired solution is a de Sitter type solution, while in region z > A it is an 
approximation to the Schwartzschield solution. Denote the constants that describe the 
desired solution in both the regions as 

a = 0; B ; b = z < A 

a; B = 0; b z>A ^ 

The lacing conditions follow from Table 1 and reduce to 




e^'J A = [e1 A = [eV]| A = 0, (9) 

that is to ensuring continuity on the discontinuity surfaces of functions 

7, 7', (3, (3' . Under assumptions (8) made, conditions (9) determine four of seven 

problem parameters p ; A ; B ; A; p; b; a. 

In what follows it would be convenient to introduce parameter zq = — a/p in lieu of 
a and use not parameters p ; A Q ; B Q ; A; p; b; z themselves, but their dimensionless 
combinations with retaining the old notations for the dimensionless parameters (excluding 
zo), 

p ^p A; p^pA; A -> b -> ^; e -> ^. (10) 

Upon the transition to the dimensionless parameters the number of the parameters 
becomes six. Thus, the problem of lacing possesses the self-similarity. 
Assume that constants p and A pl are 

po = 2.77, A p 2 = 6.94376 . (11) 
From conditions (9) we find: 

p = 2.7887, e = 0.954694, B = 2.8803 • 10~ 4 , b = -0.305319. (12) 

Hence, the accurate solution to the centrally symmetric static problem with the dis- 
continuity surface has been constructed, such that the solution is an analog of the de 
Sitter solution on the left of the surface and of the Schwartzschield solution on the right. 

Energy-momentum tensor (2) can be analyzed by the well-known Eckart method. In 
gauge 

A = Const, A^ = 0, (13) 

which is valid in the centrally symmetric static problem as well, the current density vector 
is written as 

j a = -2\A a . (14) 

In the geometrodynamics scheme under discussion, the retained substance density p 
is related to vector j a as 

r = P -u a . (15) 

In the centrally symmetric problem, 
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u a = (exp (-7/2) ; 0; 0; 0) ; u a = (-exp ( 7 /2) ; 0; 0; 0) . 

The determined energy density U and pressure P result in the following expressions 
for these quantities: 

-U = T ° = 3exp (-7) • 2 + A. (16) 

P = i (T^ + Tl + T|) = e:rp (-7) • 2 + A. (17) 

From (16), (17) it follows that for the static centrally symmetric problem the equation 
of state of geometrodynamic medium is 

p = (-5/9) U. (18) 

Assume that for the geometrodynamic medium the law of degradation of energy is 
valid in the form, in which it takes place for viscous heat-conducting continuum. This 
assumption leads to Maxwell cross relation 

d_ (dS\ = d_ (d_S_\ 

dV\dP) dP\dV)' 1 ; 

Having solved the differential relation appearing from (19), we arrive at the relation 
of temperature T to specific volume: 



(here Qo is some dimensional constant), and entropy density S is 

32/3 

S = ^A +S °' (21) 
In the stationary centrally symmetric problem under discussion in the region gauge 
(13) is used in all the thermodynamic quantities are constant and expressed in terms of 
A. 

The isentropic sound speed for the internal part of the solution is 

c 2 s = -4=. (22) 



From the condition of equal sound speed and light speed we obtain relation 

c\ = 87T70. (23) 

From (22), (23) it follows that the criterion value A cr of A, on whose achievement the 
discontinuity surface appears, is 

7 2 

Acr = 27.33.7r2.rn 2 ' (24) 
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The value of A = A cr is established as a result of some dynamic transient, whose 
description is beyond the scope of this paper. The transient apparently proceeds until A 
in the internal part of the solution rises up to a level, at which the sound speed becomes 
equal to the light speed. 

The A is related to the parameters (10) as A = B Pq and the r is nothing else but 
r = —b. Relation (24) for the solution under discussion is therefore equivalent to 

B °* 2 = (25) 
Thus, the condition of equal sound speed and light speed leads to the relation between 
the two independent constants po, A pl. The dependence of A Q pl on p following from 
condition (25) is broken at p = 2.77, at higher p equation (25) has no solutions. In the 
range of the values of p$, at which equation (25) is solvable, the curve —b (po) monoton- 
ically decreases and reaches its minimum at the boundary value p = 2.77. In so doing 
A pl = 6.94376, b = —0.305319. It is these boundary values of the constants p , A pl 
that are used in (11). 

The existence of the minimum in —b (p ) means that the solution with constants 
(11) possesses stability in the following sense. If there were no minimum, the parameter 
—b (po) could decrease to zero; as the parameter is equal to the Schwartzschield radius, 
the vanishing would essentially mean the mass "vanishing" . 

With gauge (13), the dynamic coefficient of viscosity is 

' = &■ (26) 

In viscous medium, the solution evolves according to the dynamic equations until 
Reynolds number R reaches a critical value Ro leading to appearance of instability of 
some motion modes, that is turbulence. 

Pose the problem to find size L , beginning with which the Reynolds number reaches 
the Rq. From the definition of the Reynolds number it follows that with the system size 
Lq, characteristic medium velocity u, and mass m observed from the outside the following 
relation takes place: 

L ° = u-l^GXm- (27) 
To find L , as it follows from (27), some assumptions of the value of u as well as an 
assumption allowing A to be estimated have to be made. 

As for u, assume that it characterizes the perturbation velocity, that is take 

u = c. (28) 
In view of this assumption, relation (27) becomes 

= (29 > 

where ro is the gravitational radius of the object. Assume that the coefficient in brackets 
is on the order of one. Then from (29) it follows that 

L ~ t~~ • (30) 
Ar 
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Find the characteristic size L for the above-constructed discontinuous stationary cen- 
trally symmetric geometrodynamic solution. We obtain 

o7.o3._2 

£o~ % 2 ° ^120-r . (31) 

Thus, the intolerance of the centrally symmetric solution to perturbations can appear 
at the values of the radial variable L much higher than the gravitational radius. However, 
its specific value depends on the choice of the Reynolds number for the criterion value. 

Hence, equations (1) with tensor T a p of form (2) admit the existence of discontinuous 
solutions. The possibility in itself to construct the discontinuity-surface solution is a 
nontrivial fact. Note that, for example, the geometrodynamic equations without A term 
do not admit construction of the solution obtained in the paper. 

The discontinuous solution constructed is regular in the entire radial variable range, 
which the geodesic completeness takes place in. What makes the solution to be rearranged 
with the transition of the Schwartzschield branch to the de Sitter branch at certain radial 
variable? The question touches upon the fundamental problem of the unlimited energy 
cumulation instability, that is the problem addressed in many papers. We do not have 
a full answer to the question, at least, for the reason that in this paper we restrict our 
consideration to the scope of the static centrally symmetric problem. However, we can 
specify the peculiar geometrodynamics features that can lead to the singularity removal 
by the discontinuity surface formation. 

First, in the scheme under discussion the space-time evolution proceeds in accordance 
with the dynamic equations at any time. No manual introduction of equations of state or 
coefficients of viscosity is admitted. As a result, meeting the energy dominance condition, 
which collapses are typically associated with, is not guaranteed whatsoever. 

Second, with increasing A viscosity decreases, Reynolds number increases. When the 
Reynolds number has reached some criterion value Ro, the solution becomes intolerant to 
the perturbations of characteristic size (31). It is hard to tell as applied to geometrody- 
namic medium at what specific value of Ro the perturbations will rearrange the solution. 
But that the solution rearrangement is inevitable follows from the fact that smaller and 
smaller-scale perturbations become unstable as the singularity is approached. At some 
stage the turbulization will come to be caused by quantum fluctuations of vacuum. 

The above-mentioned features of the dynamic equations allow us to suggest that the 
mechanism of the solution regularization in geometrodynamics through the discontinuity 
surface formation can be of general nature. This paper has validated this suggestion by 
the specific example. 

The work was carried out under support in part by the International Science and 
Technology Center (Project # KR-677). 
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